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Abstract 

On a cubic lattice the zero-momentum meson states have one of 20 possible A^*^ combinations 
where A labels the irreducible representation of the octahedral group. Each continuum bottomo- 
nium state with specific J^^ quantum numbers is contained within one or more of the lattice A^*-" 
states. In this work, bottomonium quark-antiquark operators are constructed for all 20 lattice 
A^*-^ combinations which allows many continuum high angular momentum states to be accessible 
as ground states of their associated lattice channels. From a dynamical simulation, realistic results 
are obtained for S-, P-, D-, F- and G-wave bottomonium states. 
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I. MOTIVATION 



Although the hghtest spin-triplet S-wave and P-wave bottomonia have been known 
experimentally for many yearsflj, the spin singlets have just been discovered in recent 
experiments 2- sl- Beyond the P waves only the T{1^D2) has been observed so far 
Meanwhile, additional resonances, which may have an interpretation as four-quark combi- 
nations, are being found with masses in the bottomonium region js, 9|. Also, observation of 
the rib{2S) has recently been claimed, though with a surprisingly small mass 10|, 111] - 

Dynamical lattice QCD computations of the masses of S-wave and P-wave bottomonia, 
including radial excitations, appear in 12|-[l6|. Although some of these publications also 



commented on D-wave states, the first lattice study of the full set of D-wave quantum 
numbers is very recent \v\ . 

There are many important reasons for extending lattice studies to higher angular mo- 
mentum states. D waves and F waves provide an opportunity to confirm and refine lat- 
tice methods, such as lattice NRQCD that uses a nonrelativistic expansion for the valence 



quarks 



18| . Quantitative predictions for the masses can subsequently be tested experimen- 



tally, and comparison of reliable lattice results with today's phenomenological models can 
offer insights into various model assumptions. In particular, the spectrum as a function of 
orbital angular momentum compared to radial excitations depends on the spatial behavior 
of the quark-antiquark interaction. Extending the study of bottomonium to higher partial 
waves augments the information gained by the observation of radial excited states. The 
lightest G-wave bottomonium is expected to lie just below the BB threshold with F waves 



further below and H waves above 



19| . so these bottomonium channels offer an opportunity 



for using lattice QCD to study threshold issues in detail, and to confirm which side of thresh- 
old the G waves are on. Furthermore, knowledge of the bottomonium masses will be crucial 
for understanding the four-quark candidates and any other exotic hadrons that could be 
nearby. Eventually, radiative transitions among these states could be considered |20i]. 

In general, a lattice operator will couple to all hadrons with the corresponding quantum 
numbers, but the lightest hadron will have the cleanest signal. Often the lightest hadron 
provides the only usable signal in a lattice study. As a result of cubic symmetry, there are 
just 20 mutually orthogonal A"^*^ options available to lattice simulations. Specifically, the 20 
options originate from the reduction of continuum angular momentum J to the five lattice 
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octahedral representations A, each with P = ±1 and C = ±1. All S-wave and P-wave 
bottomonium states are included in that list of 20 "ground states", along with just 10 higher 
waves. All other continuum bottomonium quantum numbers are only accessible as excited 
states in the 20 channels. As well, any hybrids, exotics, and multiquark states would also 
be included within the list of 20 possible lattice channels. 

It is not known, a priori, which hadron is actually the lightest in a given A^*^ channel. 
For example, the quark model suggests that a 9+^ bottomonium state will be the lightest 
in one of the 20 channels, but that same channel will couple to a 0^^ exotic hadron if one 
exists. To discover which is the lighter state, it is necessary to perform an explicit lattice 
simulation and analysis. Identifying the lightest mass in each of the 20 channels is a natural 
starting point for lattice QCD explorations of bottomonia at higher angular momentum. 

The specific contributions of the present work are twofold. First, a set of quark- antiquark 
creation operators for all 20 combinations of lattice quantum numbers is constructed and 
the lowest continuum angular momentum is identified. Secondly, numerical simulations are 
used to show that many of the higher angular momentum bottomonium states are accessible 
using present lattice QCD methods and existing gauge field configurations. 



II. BOTTOMONIUM OPERATORS FOR ALL LATTICE QUANTUM NUMBERS 

Conservation of angular momentum is a consequence of invariance under spatial rotations, 
and is realized as a tower of representations labeled by non-negative integers J G {0, 1,2,.. .}. 
On a spatially cubic lattice, this symmetry is reduced to the five irreducible representations 
A G {Ai, A2, E, Ti, T2} of the octahedral groupj21|. The well-known correspondence between 
J and A is given in Table [H In principle a sophisticated analysis of excited states in lattice 
data for each A channel would allow the unambiguous identification of J for each observed 
lattice state. However, experience from basic quantum mechanics and the quark model 
suggests that larger J leads to larger bottomonium mass, so in practice it is simple to infer 
each J value from a single channel (unless a non quark model hadron appears with a nearby 
mass) . 

Total angular momentum is the sum of orbital and spin components. The orbital an- 



3 



Table I: The number of copies of each octahedral irreducible representation A, of dimension d, for 
the subduced continuum representation J. 

A{d) J 

0123456789--- 

^i(l) lOOOlOlOll--- 

^2(1) OOOlOOllOl--- 

E{2) 0010111121--- 

Ti(3) 0101121223--- 

T2{3) 0011112222--- 

gular momentum L of each operator corresponds to the number of uncontracted^ spatial 
derivatives, but it is important to remember that A (the lattice version of J not L) is the 
conserved quantum number. An operator intended for a specific L will couple to other L 
values too, unless the fields are smeared in a way that purposefully eliminates that mixing. 

One operator for each of the 20 A^*" channels is provided in Table [Tll They are written in 
terms of Pauli matrices (Xj and the symmetric lattice covariant derivative Aj which appears 



beyond the P waves in the following forms: 

Dij = ^(A,A, + A, A,) - hij{Al + AI + Aj) , (1) 

Dijk = ^(AjAjAfc + permutations) , (2) 

Dijki = ^(AjAjAfcA; + permutations) , (3) 

Dijkim = ^(AjAjAfcA^Am + permutations) , (4) 

Dijkimn = ^(AjAjAfcA,A„A„ + permutations) , (5) 

Dijkimnopq = ^ ( Aj Aj Afc A; A„ A„ A^ Ap Ag + permutatious) . (6) 



The entries in Table [IT] can be verified by using standard group theory methods to calculate 
multiplicities from the character of each conjugacy class, but an intuitive understanding 
can be obtained more readily. For example, the operator is listed as D123. Any spatial 



A lattice Laplacian, for example, does contain a Lorentz contraction and does not alter the orbital angular 
momentum. 
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rotation of this operator on a lattice will reproduce the same operator except that its indices 
may be permuted and for certain rotations there will be an overall minus sign due to the 
vector nature of a spatial derivative. These spatial rotations have the properties of the 
one-dimensional representation A2, including the overall minus signs for the appropriate 
rotations. P and C properties can also be verified. 

For use with lattice NRQCD, each of the operators in Table [ITl will be sandwiched between 
a bottom creation operator (x) and an antibottom creation operator xK^)^ then summed 
spatially to produce a zero-momentum interpolating operator: 

J2^\x)nxHs)- (7) 

X 

The ip'^i^x) and xK^) '^^^ contain internal summations as well, to represent spatial smearing. 
In our simulations, covariant derivatives act on the field to the right and not to the left. 
Charge conjugation is maintained exactly for all operators because of the sum over x in 
Eq. ©. _ 



A table similar to Table [ITl appears in [22], but with a few notable differences. In 22 1 
the authors listed only S, P and D waves, and their double derivative operator Ajj is not 
symmetrized in its indices, whereas we prefer to use a symmetric definition so exact sym- 
metries are apparent. For the E^'^ channels, which are two-dimensional representations, 
the authors of 22| chose to list three nonindependent projections whereas we prefer to list 
the two independent projections. They also included three additional operators in the 
and T2 channels, intended to probe D-wave states which, however, are not expected to be 
ground states in their respective A^*^ channels; these operators are displayed in our Table UTT] 
for convenient reference. For relativistic operators rather than NRQCD operators, as might 
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24|. 



be more appropriate for charmonium, see 

Although the list of 20 quantum numbers A^*^ is complete, the operators listed for each 
are not unique. The operators displayed in Table [ITl were chosen with standard bottomonium 
states in mind, so each contains one quark, one antiquark and L derivatives. These operators 
will in general have some overlap with exotics, hybrids and multiquark states as well, but a 



search meant to target those states should construct other operators with t 



le goal of having 



a larger overlap with the particular states of interest (see, for example. 

The ^Di meson does not appear in Table HT] because it has the same J^'-' as the lighter ^5*1 
meson. Similarly, the meson does not appear in Table [ITl because it has the same 

jPC 

as 



Table II: Creation operators 0, for each irreducible representation A of the octahedral group with 
specified parity P and charge conjugation C. Columns labeled J^*^ and "^^^^Lj contain the quan- 
tum numbers of quark-antiquark states obtained using the standard relation of S and L with parity 
and charge conjugation. 



^PC jPC 


25+1^^ 


n 


^1 


u 




1 


T — 


1 — 




{cTi, (72, Cs} 




1+- 


'Px 


{Ai, A2, A3} 


^1 


n++ 




AifTi + A2cr2 + AafJs 


^1 


1++ 


^Pl 


{A2Cr3 — A3(T2, A3C71 — AiCr3, Ai(T2 — A2(Tl} 


11/ 


+ + 
L, 


'^2 


{(AiCJi - A2Cr2)/\/2, (AifJi + A2Cr2 - 2A3a3)/\/6| 

_LJ_ £^ £^ / I y 1 \ J__L' ^ ^J^J//V J 




+ + 


^P2 


{A2Cr3 + A30"2, A30"l + Ai(T3, Ai(T2 + A2Cri} 


11/ 






{(Dll - D22)lV2, (Du + Z)22 - 2L»33)/^/6| 








{1)23,^31,^12} 


TP 


9 




{p230-l - Di30-2)/\/2, (-D23O-I + 1)31(79 - 20x20^)1 


T — 


2 — 




1(^22 — -D33)cJl + Dl3Cr3 — I)l2C"2, (D33 — i?ll)o"2 + -D2I 0"l — -D23C3, 
{Dii - 022)^3 + D^2<y2 - ^3lO"l} 


^2~ 


3 




^12C"3 + 023(71 + -D3l(72 




3+- 


^3 


-D123 


n- 


3+- 




{1)122 — -C>i33, D233 — D2ll,D3u — D322} 




3++ 


'Fs 


{D221 - D33i)ai + (-D332 - -Dll2)f72 + (^113 - ^223)0"3 




4-+ 


'G, 


{-^2223 — -^3332,-03331 — -Diii3,Diii2 — -D2221} 




4— 


'G, 


(-C2223 — -C3332)o"l + (-D3331 — -Diii3)(J2 + (-D1112 — -02221)173 


E+- 


5+- 




{(-C23III - -Cl3222)/\/2, (-D23III + -D13222 - 2.Di2333)/\/6} 




6-+ 




-D112222 + -D223333 + -D331111 — -D221111 — -D332222 — -D113333 




9+- 




-C122233333 + -C2333IIIII + -D311122222 — -C133322222 — -D211133333 — -C3222IIIII 



the lighter ^P2 meson. Neither of those omissions is a consequence of lattice discretization, 
since even the continuum quantum numbers are shared. On the other hand, the absence of 
the meson from Table [ITl is a consequence of lattice discretization: its continuum quantum 
numbers (4"''+) are not shared by any lighter state, but its lattice quantum numbers {Af'^, 
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Table III: Creation operators that duplicate some quantum numbers from Table |TT1 but could be 
useful for studies of the D-wave spectrum. They were first listed in [22 1. 



^PC jPC 2S+lLj ^ 

1 '^Di {Diiai + Di2(J2 + Di3a3,D2icri + D22cr2 + -C»23cr3, D^iai + 032(^2 + -Dsso-s} 
3— 3/^3 {SDnai - 20^2^2 - ^Disa^, 3D22CT2 - 21)23(^3 - 2D2iai, 
- 2D3iai - 2032(72} 

T^ - 3— 32^3 {{D22 - D33)(Ji + 2Di2(J2 - 2D13CJ3, {D33 - Dn)(T2 + 202303 - 2D2iai, 

{Du - D22)(J3 + 21)31^1 - 2Z)32fT2} 



E^^ , Ti^, T2^) are all shared by the ^Pj mesons. Therefore lattice signals of the ^-Di, 
3^2, '^^4 and many similar cases must be sought as excited states in the corresponding A'^'" 
channels, unless operators can be tuned carefully enough to remove any coupling to the 
lighter mesons in each channel. Analysis of excited states in any A^*" channel is beyond 
the scope of the present work with one exception: we will address the T]h{2S) because of the 
present experimental situation 1^, 11 1. 



III. LATTICE SIMULATIONS 



All numerical simulations performed for this work utilize a 2+1-flavor dynamical gauge 
field ensemble generated by the PACS-CS Collaboration[25] and made available through 
the Japan Lattice Data Grid 26j as part of the International Lattice Data Grid 27]. This 
ensemble uses the Iwasaki gauge action and the nonperturbatively tuned clover action for 
up, down, and strange quarks. The light quark (up and down) masses are such that the 
pion mass is 156 MeV/c^, i.e. near the physical value. There are 198 configurations with 
lattice volume 32^ x 64. The four input parameters, along with four basic output quantities, 
are displayed in Table llVi Not provided by the PACS-CS Collaboration but needed for the 
present work is the mean link in Landau gauge that was computed 20| to he ul = 0.8463. 

The bottom quark and antiquark are described by the lattice NRQCD Hamiltonian jl8 1 
omitting terms at 0{v^) and beyond. Specifically, the Hamiltonian from the Appendix of 



28| is used with Cj = 1 for z < 6 and Cj = for z > 7, with tadpole factor Uq = ul = 0.8463 



and with stability parameter n = 4. 



7 



Table IV: Input parameters for the ensemble of gauge configurations generated by the PACS-CS 
Collaboration[25]. Some standard output quantities provided by PACS-CS are also listed. 



Inputs Standard outputs 

f3 = 1.90, a = 0.0907(13) fm 
Kud = 0.13781, = 156(7) MeV/c^ 
Ks = 0.13640, rriK = 554(8) MeV/c^ 
csw = 1-715, (□) = 0.570870(5) 



In previous work[2o| the bottom quark bare mass was tuned to match the experi- 
mental T]b{lS) mass, with result Mb = 1.945. However, a recent measurement by the 
Belle Collaboration[29] is significantly larger than the rii,{lS) mass reported by previous 
experiments |l|. Tuning to the Belle result would require a slight increase in the bottom 
quark bare mass of lattice NRQCD. To avoid this unresolved issue, it is preferable to use 
the well-established T(IS') mass|l| for tuning instead. This also leads to a slight increase in 



20| . and for the present project we 



the bare mass parameter relative to our previous work 
use Mb = 1.95. 

To reduce statistical fluctuations, it would be optimal to average over results obtained 
from putting the creation operator at each lattice site. A practical implementation of this 
general idea is employed here: heavy quark propagators are constructed from a random-U(l) 



wall sourcej30| spanning all spatial lattice sites at a flxed Euclidean time. A further reduction 
of statistical fluctuations can be achieved by selecting another Euclidean time on the same 
gauge conflguration, separately inserting a random-U(l) wall source there, and averaging the 
two correlation functions on that conflguration. For the present work, averages are taken 
over 64 separate calculations of all correlation functions per conflguration, corresponding to 
a random-U(l) wall source at each Euclidean time on the 32^ x 64 lattice. 

Excellent ground state signals are obtained for S and P waves with a single correlation 
function using a local operator. As one goes to higher waves, simulation energies increase 
and correlation functions decrease more rapidly and become noisier due to the presence 
of more derivatives in the operators. It is advantageous to utilize more correlators in each 
channel and to smear the operators to enhance the ground state contribution at smaller time 
separations. Beyond the P waves, gauge-invariant smearing |i31j] is applied to the bottom 
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quark operator, 

^(x) ^ (1 + 0.15A2)^Xx) . (8) 

Three options are considered: s = (local), s = 1 (smeared), and s = 2 (doubly smeared). 
In every case the antiquark field x{x) remains unsmeared. 

With local, smeared and doubly smeared operators available for both creation (s) and 
annihilation (s'), each A^*-^ has nine correlation functions to be analyzed simultaneously. 
The nominal fit function is 

N 

9s,At -to) = Y. fs{n)fAn)e-''-^'-'-^ (9) 

n=l 

where the number of exponentials is increased until the x^/d.o.f. is unity or less. En 
denotes the simulation energy for the nth meson while fs{n) and fs'{n) denote the coupling 
of the nth meson to the creation and annihilation operators respectively. 

In practice, the lattice data cannot always resolve individual higher-energy states and that 
requires a generalization of Eq. dH]). As a simple example consider two operators, A and B, 
that each couple to two physical mesons having energies Ei and i?2- Suppose these energies 
are unequal but not distinguishable with the numerical precision of the available lattice data, 
so they get represented by a single exponential during the fitting process, Ei ^ E2 ^ -Efit- 
The nominal fit function (with to = 0) would give 

{0\AA^\0) = [fAil)fe-^^' + [fAi2)fe-^^' ^ ( [/^(l)]^ + [/a(2)]^ je"^-*, (10) 
{0\BB^O) = [/B(l)]^-^^* + [/i,(2)]^-^^* ^ ([/B(l)]' + [/B(2)]')e-^-*, (11) 
{0\BA^O) = ^(l)/^(l)e-^^* + /^(2)/^(2)e-^^* ^ (/^(1)/b(1) + /^(2)/^(2))e-^-*.(12) 

The three coefficients of e~^^* contain only two degrees of freedom as assumed in Eq. (|9]). 
The same is true for the coefficients of e"^^*. However, the three coefficients of e~^"'* do not 
contain just two degrees of freedom. To account for this potential inability of lattice data 
to resolve excited states, Eq. (Q can be generalized to 

AT' AT 

<7.,,,(t-to) = J]/.(n)/,,(n)e-^"(*-*«)+ /^y (n)e-^"(*-*°) (13) 

n=l n=N'+l 

which says the lowest A^' states are resolved individually but higher states may not be. 
Confident extraction of the ground-state energy requires A^' > 1, which is satisfied easily by 
all of the data discussed in this paper. 
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The methods described above are sufficient to permit explorations in the S, P, D, F, 
G, H and I channels as presented in the following section, with reasonably precise results 
obtained in the S, P and D channels. Since F and G waves have not been discussed in 
p..evi„„. lattice s„„u.ati„„s but a.e expected to Me below the BB thtesholdQ. we focus 
particular attention on the F- and G-wave channels by implementing, in addition to the 
methods described above, stout spatial links U as proposed in 32|: 



U f/(^) ^ f/(2) ^ y = U , (14) 

f/("+i) (x) = exp (^g(") (x)) f/(") (x) , (15) 

where Q^^\x) contains a pro duct of gauge links weighted by a parameter p. [For the precise 



definition of Q^f!^\x), see |32|.] Our F- and G-wave simulations use p = 0.2 and Up = 16 for 
source and sink operators, but no stouting within the NRQCD propagator. 

Lattice simulations involve both systematic and statistical errors. In the present context, 
systematic errors arise from the neglect of 0{v^) terms in the NRQCD Hamiltonian, the 
neglect of 0{as) corrections to the NRQCD coefficients q, the nonzero lattice spacing, the 
finite lattice volume, and up and down quark masses that are slightly larger than their phys- 
ical values. Although we do not have any quantitative means to determine systematic errors 



for ourselves, the authors of 12| studied bottomonium with the same truncation of NRQCD 
terms and tadpole- improved tree-level Cj definitions on lattices having discretizations and 
volumes similar to ours, though their up and down quarks were significantly heavier than 
ours. The thorough discussion in [12] concluded that systematic errors in typical T mass 
splittings are at the few-percent level. We take this as a generic expectation for our sys- 
tematic errors too, though only statistical errors will be reported below. Interestingly, the 
statistical errors will turn out to be at the few-percent level also. 



IV. LATTICE RESULTS 



To provide a visual impression of the quality of lattice data, sample D-wave and F-wave 
correlation functions are displayed in Figs. [1] and [2] respectively. Numerical results for the 
ffis are listed in Table IVl Because the smearing parameters were not optimized individually 
for each operator, two of the ffis do not beneffi from using all nine combinations of source 
and sink smearings; Table |V] records the number of combinations used in each case. Fits 
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Figure 1: Correlation functions for E with equal smearing at source and sink. The lightest meson 
in this channel is ^Z)2 bottomonium. 



employ the standard correlated method 33|. With a source at Euclidean time to, the S- 
and P-wave fits include 1 <t — Iq < 31, D-wave fits use 1 < t — Iq < 27, and all higher- wave 
fits use 1 < t — to ^ 23. The conversion from the lattice simulation energy i?sim to a meson 
mass in physical units (also shown in Table |V]) is accomplished by recalling that Mf, was 
tuned to the physical T(IS') mass. Explicitly, 



Mc' = M(T(15)e.periment)c' ^ ^ (^^'"^ ' ^sim(T(15)) 



(16) 

a \ / 
where the lattice spacing, a, is taken from Table IIVI 

The T]h{lS) mass in Table IVl is consistent with the Belle measurement and larger than 
the PDG average: 



Miri,{lS)) 



(17) 



9390.9 ± 2.8 MeV/c2, PDG average[lj 
9401.0 ± 1.9^^:^ MeV/c^, Belle[29] 

9404.0 ± 0.8 MeV/c^, this work, statistical errors only. 
It should be noted that this lattice result is consistent with some lattice determinations that 



predated the Belle measurement 
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20|. Studies of specific aspects of this issue are also 
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Figure 2: Correlation functions for T2 with equal smearing at source sink. Stout links are used 
at source and sink as well. The lightest meson in this channel is bottomonium. 

available in QQ- 

The P-wave and D-wave states are in reasonable agreement with the available experi- 

n 

mental measurements, but splittings among D waves are not resolved, in contrast to |17| . 
The three F-wave results from Table |Vl 



MCFs) 



10.369 ±0.017 GeV/c2 
10.383 ± 0.016 GeV/c^ 
10.367 ±0.018 GeV/c^ 



(18) 
(19) 



are in excellent agreement with a model prediction of 10.37 GeV/c^[19]. The two G-wave 
results from Table |Vl 



MCd) = 10.581 ±0.017 GeV/c^ 
= 10.587 ± 0.018 GeV/c^ 



(20) 
(21) 



are about one standard deviation above the experimental BB threshold at 10.56 GeV/c^. 
Recalling that the quoted lattice uncertainties are statistical only, our results are also rea- 
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Table V: Simulation energies and resulting bottomonium masses obtained from fits to lattice data. 
Also shown are the number of correlators (i.e. smearing combinations) used in the fits and the 
number of exponentials needed to reduce the x^/d.o.f. to unity or less. All operators are from 
Table HTl except the one labeled ^D^) which is from Table UTTl Errors are statistical only. 



j^PC 








mass [GeV/c^] 


^1 


1 


4 


0.23420(6) 


9.4040(8) 


T-, — 
1 


1 


4 


0.25993(7) 


9.460 


T+- 
^ 1 


1 


4 


0.451(4) 


9.876(11) 


1 


1 


4 


0.438(4) 


9.847(11) 


T++ 
1 


1 


4 


0.449(4) 


9.871(11) 


E++ 


1 


4 


0.457(7) 


9.889(17) 


T++ 

2 


1 


4 


0.454(7) 


9.882(17) 


E-+ 


9 


5 


0.591(5) 


10.180(16) 


2 


9 


7 


0.599(4) 


10.198(14) 


E — 


9 


6 


599(4) 


10 lQ8fl4) 


Z 


9 


6 


0.588(3) 


10.174(13) 




9 


6 


0.600(3) 


10.200(13) 




9 


6 


0.601(3) 


10.202(13) 


rn+- 


9 


6 


0.684(4) 


10.369(17) 


^2 


9 


5 


0.690(3) 


10.383(16) 




9 


5 


0.683(5) 


10.367(18) 




6 


4 


0.775(1) 


10.581(17) 




6 


4 


0.778(3) 


10.587(18) 



sonably close to the model prediction of 10.52 GeV/c^ 19| (which, however, lies below the 
BB threshold). 

Fits to lattice data in H-wave and I-wave channels could only support three or four 
exponentials, had no stout links, and lay above the BB threshold, so meaningful mass 
predictions were not obtained. We did not attempt an L-wave simulation. 

The behavior of orbital excitations compared to radial excitations depends on the form 
of the potential. Figure [3] shows the qq potential used in the quark model of Godfrey 
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Figure 3: Main graph: The qq potential for mesons in the relativized quark model of Godfrey and 
IsgurjsSl; see their Fig. 12. Upper inset: The spectrum of low- lying energy levels from a simple 
Coulomb potential. Lower inset: The spectrum of low- lying energy levels from a simple linear 
potential. 

and IsgurjsSi, along with their determination of the location of some T excitations in that 
potential. The T(IS') is in the Coulomb regime but higher radial excitations move into 
the transition region between Coulomb and linear. The insets in Fig. [3] show radial and 
orbital excitations in the simple cases of purely Coulomb and purely linear potentials. In 
calculations and from experimental data one finds that the ground-state P-wave states have 
masses not very different from the first S-wave radial excitations. There is at least some 
resemblance to the expectation from a Coulomb-like potential. The situation changes as 
angular momentum is increased. At G waves the ground state would be comparable to 
the T(55') in the Coulomb limit but comparable to T(35') in the linear limit. The lattice 
spectrum from the present work is displayed in Fig. IH where the G wave is seen to be 
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Figure 4: Lattice results for the bottomonium masses (data points with statistical error bars) 
compared with experimental data|l| (horizontal lines). Radially excited T masses are shown for 
comparison to orbital excitations. 

comparable to the T(4S') mass, i.e. intermediate to the Coulomb and linear cases. For the 
lattice results, the G — F splitting is comparable to the F — D splitting and thus more 
reminiscent of a linear-dominated potential than a Coulomb-dominated potential. Perhaps 
the most important conclusion about G waves in this work is that statistical errors are small 
enough to warrant future systematic lattice study to determine, for example, whether the 
G waves lie above or below threshold. 

Although this work has focused exclusively on the ground state masses for each A^*^, 
excited state masses were obtained during the fitting process. In light of the present ex- 



perimental situation 



10 



the r]i,{2S) is of particular interest. Our present work relies 



on local operators for the S waves, but our previous workj20j included smearing and thus 
produced a more precise result for the rih(2S) mass. Since it was not reported explicitly in 
that publication, we include it here: 

M (T(25)) - M {r]b{2S)) = 24 ± 3 MeV/c^ (22) 

where the error is statistical only. This hyperfine splitting is significantly smaller than the 
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recently announced experimental measurement 10|,llli, M {T{2S)] 
MeV/c^, but our result is consistent with other lattice predictions 

3|. 



-Mjj]h{2S)) = 48.7±2.7 
I2I . Il4j-|l6| as tabulated in 



V. CONCLUSIONS 

The study of bottomonium at increasing angular momentum provides knowledge that 
augments that obtained from the spectrum of radial excitations. In this paper we provide a 
set of operators in lattice NRQCD for higher angular momentum bottomonium states and 
explore the feasibility of their use in numerical simulations. 

A set of quark- ant iquark operators for all lattice irreducible representations, A^*^, has 
been given in Table [Til In principle, this set of operators allows for bottomonium states 
with 16 distinct continuum J^'-' combinations to be described as ground states in their 
respective lattice channels. All the S- and P-wave states are in this set. However, for 
higher angular momentum, some continuum J^^ combinations would only appear in 
channels where a lower angular momentum state is present. These states would have to be 
extracted from lattice simulations as excited states. Since ground state energies are much 
easier to determine from lattice simulations, the states described by the operators in Table HTl 
are a natural starting point for numerical studies of the high angular momentum states of 
bottomonium. 

Our numerical simulations produced clear signals for S-, P-, D-, F- and G-wave states 
as displayed in Fig. HI These are the first lattice simulations of F-wave and G-wave bot- 
tomonium, and they are near quark model expectations. The G waves are very close to 
the physical BB threshold, and statistical errors in the lattice simulation are small enough 
to encourage future studies of the systematic errors so that lattice methods can ultimately 
predict with confidence whether the G waves lie above or below threshold. Table HTl contains 
ground-state operators for H, I and L waves also, though these will be significantly above 
the BB threshold. 
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